Virtual Photon Correction to the K + — » 7r + 7r°7r° Decay 
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Abstract 

We consider electromagnetic corrections to the non-leptonic kaon decay, K + — ► 7r + 7r 7r , due to 
explicit virtual photons only. The decay amplitude is calculated at one-loop level in the framework 



of Chiral Perturbation Theory. The interest in this process is twofold: It is actually measured by 
the NA48 collaboration from one side, and, the value of the amplitude at the 7T7T threshold gives 
access to 7T7t scattering lengths from the other side. We found that the present correction is about 

T 1 

( JC ■ 5 to 6% the value of the Born amplitude squared. Combined with another piece published recently, 

^ this fixes the size of isospin breaking correction to the amplitude squared to 7% its one-loop level 

^o : 

(^i value in the absence of isospin breaking and at the center of Dalitz plot. 

o 
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I. INTRODUCTION 

Starting May, Nicola Cabibbo submitted to the arxiv a letter proposing a "potentially 
accurate" method for the determination of tttt scattering length a{] — Oq fr° m the 7T°7r 
spectrum in K + — > 7r + 7r°7r° near the tt + tt~ threshold |l|. Next day, Hans Bijnens and Fredrik 
Borg submitted a paper treating isospin breaking in K — ► 3tt decays in the framework of 
chiral perturbation theory [2|]. In their study, the authors of Ref. |2] took into account 
strong isospin breaking and local electromagnetic corrections. The former is proportional 
to the up and down quark mass difference m u — m&. With respect to local electromagnetic 
correction, it includes the effects of mass square difference between charged and neutral 
pions M^ ± — M^ and electroweak counterterms. 

Let us justify the necessity of the present work. The NA48 collaboration investigates 
actually the CP violating asymmetry in K^ — > -k^-k ^ decays. So, there is no need in 
principle to upgrade the experimental setup in order to study the 7r°7r° spectrum near the 
tt + tt~ threshold. From the theoretical point of view, the tttt scattering length Oq — clq is 
very sensitive to the value of the quark condensate. Hence, an accurate measurement of 
this quantity would allow a best understanding of the quantum chromodynamics vacuum 
structure. 

The most convenient theoretical framework to study non-leptonic weak kaon decays is 
chiral perturbation theory extended to include non-leptonic weak interactions J3| . The effec- 
tive Lagrangian constructed in Ref. J2| has been used in J4j,lijla] to calculate K ^ 3tt decays 
to next-to-leading chiral order. The theoretical prediction did not match the experimental 
data. Naturally, the need to push the chiral expansion up to higher orders is imperative. In 
a parallel direction, one has to improve the theoretical prediction by evaluating the size of 
isospin breaking corrections. To do so, electromagnetism should be incorporated appropri- 

n 

ately in the effective Lagrangian formalism [7(. 

We shall distinguish between two types of virtual photons: hard and soft. Hard pho- 
tons are integrated out from the effective theory leaving local electromagnetic interactions 
weighted by order parameters. The latter are counterterms in the sense that they absorb 
ultraviolet divergences produced by soft virtual photon loops. The isospin breaking correc- 
tions calculated in Ref. J2| are due to m u — rrid and to hard virtual photons. Obviously, they 
are ultraviolet divergent. Moreover, their size is about 10% the value of the one- loop level 



squared amplitude in the absence of isospin breaking. The authors of Ref. |2j concluded 
that the correction they obtained is not "quite enough to solve the discrepancies" between 
theory and experiment and that the effects of soft virtual photons should be included. The 
aim of the present work is to evaluate these effects. 

II. KINEMATICS 

The decay process 

K+(k) - vrVKVK+fe) (1) 

is studied in terms of three scalars 

Si = {k-pif, i = 1,2, 3, (2) 

satisfying the on-shell condition 

J2 S i = M K± + M n± + 2M n° ■ (3) 

% 

The differential decay rate is defined by 

dr = «b |A1|2d *' (4) 

with the differential phase space 

d 3 Pi 



d$ = (2tt) Wj> - fc J JJ 



n o ' 



(27r) 3 2p! 
and the decay amplitude M.. The spectrum of the phase space with respect to s 3 reads 



(5) 



d$ J_ \V 2 (s 3 ,M 2 ,M 2 ) \V 2 (s 3 ,M 2 ± ,M 2 K± ) 

ds 3 " 2V s 3 M| ± U 

where A stands for the Kallen function 

\(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . (7) 

The physical region is bounded by 

(M w± + M n o) 2 <C Sl , s 2 < {M K± - M n o) 2 , 4M 2 <: s 3 ^ {M K± - M n ±) 2 . (8) 

In what follows, we set M^± = Mk and M n ± = M^o = M n . This is safe since e 2 M 2 ± = 
e 2 Ml + 0(e A ). 



III. THE DECAY AMPLITUDE 

The decay (£Q) can proceed locally. It can also be mediated by pion or kaon multi-poles. 
Let us then write the decay amplitude at one-loop level as J9| : 

Z Z 

M = T+ Ml-Ml + {Ml-Mlf (9) 



We follow notations of Ref. [2| for the strangness-changing Lagrangian parameters. The 
decay amplitude takes then the following form: 

X = CG 8 X {8) + CG' 8 X {8,) + CG 27 X {27) , X = M,T,Z, Z. (10) 

Moreover, the subscripts and 1 will be assigned to Born amplitudes and one-loop explicit 
virtual photon corrections, respectively. 

A. Born amplitude 

The lowest order in the perturbative expansion of the amplitude comprises five diagrams 
as sketched in Fig. 1. We present the result in terms of scalar products avoiding any 
ambiguity 

To = o ( k ' Pi + k ' P2 - 4fc ■ p 3 - pi • p 3 - p 2 ■ p 3 ) , (11) 

TP = \m 2 k , (12) 

T (27) = - (2fc • Pl + 2k ■ p 2 - 8k ■ p 3 - 30 Pl • p 2 + 13pi ■ p 3 + I3p 2 ■ p 3 ) , (13) 

y 

Zf ] = -\{Ml-Ml)x 

{Ml + 2k- Pl + 2k-p 2 -4k-p 3 + 4pi • p 2 - 2 Pl ■ p 3 - 2p 2 ■ p 3 ) , (14) 

Zp = - l -{M 2 K -Ml)M 2 K , (15) 

if) = -y s [±Ml{Ml-Ml) 

+ (8M 2 K + 7Ml)(k-p 1 + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 -p 3 -p 2 -p 3 )] . (16) 

Note that Ai K vanishes if we add all contributions. 



B. Virtual photons 

In this section, we shall calculate one-loop diagrams with a virtual photon exchanged 
between two meson legs, one meson leg and a vertex, and two vertices. The various topologies 
are depicted in Figs. (2) - (9). The loop functions figuring in the expressions are tabulated 
in the appendix. A fictitious mass m 7 is attributed to the photon in order to regularize 
infrared divergent integrals. For completeness, we calculate the amplitudes proportional to 
G' 8 . Only those proportional to Gs and G 2 j will be taken into account numerically. We shall 
give the contribution of each figure separately. 

Figure 2 

We have 

Tl 8) = j{2A (M„) + 2A (M K ) 

-2{47i)- 2 {k ■p 1 + k-p 2 -4k-p 3 -p 1 -p 3 -p 2 -p 3 )[l + ln(ra 7 ) - ln{M n M K )} 
+ (8M 2 + k ■ Pl + k ■ p 2 - 4k ■ p 3 - Pl ■ p 3 - p 2 ■ p 3 )B {~p 3 , 0, M n ) 
+2{4M 2 -k-p 1 -k-p 2 + 4k-p 3 +p 1 -p 3 +p 2 - p 3 )B 1 (-p 3 , 0, M n ) 
+(8M| + k ■ pi + k ■ p 2 - 4k ■ p 3 - pi • p 3 - p 2 ■ p 3 )B (-k, 0, M K ) 
+2{4M 2 K -k-p 1 -k-p 2 + 4k-p 3 +p 1 -p 3 +p 2 - p 3 )B 1 (-k, 0, M K ) 
-(k ■ pi + k ■ p 2 - 4k ■ p 3 - pi ■ p 3 - p 2 ■ p 3 )B (p 3 - k, M w , M K ) 
+4(k ■ p 3 ) (k ■ pi + k ■ p 2 - 4k ■ p 3 

-Pi -P3-P2- P3)C (-p 3 , -k, m 7 , M w , M K )} , (17) 

TP = ^M|{-2(47r)- 2 [l + ln(m 7 )-ln(M w M i ,)] 

+B (-p 3 , 0, M w ) - 2Bi(-p 3 , 0, M w ) + B (-k, 0, M K ) - 2B x (-k, 0, M K ) 
-B (p 3 - k, M n , M K ) + 4(k ■ p 3 )C (-p 3 , -k, m 7 , M n , M K )} , (18) 
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Tf 27) = e -{4A,{M 7T )+4A {M K ) 
y 

-2(4ir)- 2 (2k -pi + 2k-p 2 -8k-p 3 

-30pi • p 2 + 13pi • p 3 + I3p 2 ■ p 3 )[l + ln(m 7 ) - \n(M n M K )] 

+ (16M 2 + 2k-pi + 2k-p 2 -8k-p 3 

-30pi • p 2 - Ylpi ■ p 3 - I7p 2 ■ p 3 )B (-p 3 , 0, M w ) 



+2(8M 2 -2k- Pl -2k-p 2 + 8k-p 3 

+30pi • p 2 - 13pi • p 3 - 13p 2 ■ p 3 )-Bi(-P3, 0, MJ 

+ (16M| - 28A; • Pl - 28k ■ p 2 - 8k ■ p 3 

-30 Pl • p 2 + 13 Pl • p 3 + 13p 2 • P3)B (-k, 0, M*) 

+2(8M| -2k- Pl -2k-p 2 + 8k-p 3 

+30 Pl • p 2 - 13 Pl • p 3 - 13p 2 • Ps)Bi(-k, 0, Ma) 

-2(fc ■p 1 + k-p 2 -4k-p 3 - 15pi -p 2 -pi-p3~p 2 - p 3 )B (p 3 - k, M n , M K ) 

-15(/c-pi + k-p 2 -p 1 -p 3 -p 2 •p3)B 1 (p 3 -k,M n ,M K ) 

+4(fc • p 3 )(2fc ■p 1 + 2k-p 2 -8k-p 3 

-30pi • j9 2 + 13pi • p 3 + 13p 2 • P3)Co(-P3, -A;, ra 7 , M ff , Ma) 

+60(/c • p 3 )( Pl -p 3 +p 2 - p 3 )Ci(-p 3 , -fe, 0, M ff , Ma-) 

+60(A; • p 3 )(fc • Pl + k ■ p 2 )C 2 {-p 3 , -k, 0, M n , M K )} . (19) 



Figure 3 
We have 



Tl s) = -T/ 27) (20) 

2e 2 
= -— {^(M,) - A (M K ) - M 2 K B {-p 3 , 0, M ff ) 

-M$B (p, M v , M K ) + M^ (A; - p 3 , M ff , M ff ) 

+(4M| + 3M 2 + k- Pl + k-p 2 -2k-p 3 

+2p x -p 2 -pi-p 3 -p 2 - p 3 )B (-k, 0, Ma) 

-(3M| + 3M* + k- Pl + k-p 2 -2k-p 3 

+2p 1 -p 2 -pi-p 3 -p 2 - p 3 )B (-k, 0, M^) 

-(M* + k-p 1 + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p 3 -p 2 - p 3 )B (p 3 - k, Mn, Ma) 

+ (M| + k-p 1 + k-p 2 -2k-p 3 

+2p 1 -p 2 -pi-p 3 -p 2 - p 3 )B (p 3 - k, M v , M v ) 

+ {Ml + Ak-p 3 )(Ml + k-p l + k-p 2 -2k- p 3 
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+2pi -p 2 -p\-p 3 -p 2 - P3)C (-P3, -k, m 7 , M„, M K ) 
+Ml(2M 2 K + 3Ml + k- Pl + k-p 2 -2k-p 3 
+2p 1 -p 2 -Pi-P3-P2- Pz)C {-k, -k, 0, M w , M K ) 
+ (M| -Ml-4k-p 3 )(M 2 K + Ml + k- Pl + k-p 2 -2k- p 3 
+2p 1 -p 2 -Pi-P3-P2- P3)C (-p 3 , -k, 0, M v , M w ) 
-Ml{Ml + k- Pl + k-p 2 -2k-p 3 
+2 Pl ■p 2 -p 1 -p 3 -p 2 - p 3 )C (0, k - p 3 , M K , M v , M w ) 
-Ml{M 2 K - Ml - 4k ■ p 3 ){Ml + k- Pl + k-p 2 -2k-p 3 
+2p 1 -p 2 -pi-p 3 -p 2 - p 3 )D (-k, -p 3 , -k, rriry, M K , M n , M w )} , (21) 

T f) = ^-M^iBoik-p^M^M^-B^M^Mx) 

-B (-k, 0, M K ) + 2B (-k, 0, M K ) - B (-p 3 , 0, M n ) 

+{Ml + k-p 1 + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p 3 -p 2 - P3)Co{-p 3 , -k, m 1 , M n , M K ) 

+ {2M 2 K + 3M 2 + k- Pl + k-p 2 -2k-p 3 

+2p x -p 2 -pi-p 3 -p 2 - p 3 )C (-k, -k, 0, M w , M K ) 

+ (M 2 K -Ml -4k- p 3 )C (-p 3 , -k, 0, M n , M n ) 

-(Ml + k- Pl + k-p 2 -2k-p 3 

+2pi -P2-P1-P3-P2- Ps)O)(0, k - p 3 , M K , M v , M w ) 

-{M 2 K -Ml-4k-p 3 ){Ml + k-p l + k-p 2 -2k- p 3 

+2 Pl -p 2 -pi-p 3 -p 2 - p 3 )D {-k, -p 3 , -k, m 7 , M K , M n , M n )} . (22) 



Figure 4 
We have 



if = -ZP = 1 Z[ 27) (23) 

= -^- M 2 K {-2A (M n ) 

+ (M| + 4Ml + k- Pl + k-p 2 -2k-p 3 
+2p x -p 2 -pi-p 3 -p 2 - P3)B (-P3, 0, M„) 



+(3M| + 2M 2 + k- Pl + k-p 2 -2k-p 3 

+2pi -P2-Pi-P3-P2- p 3 )B (-k, 0, M w ) 

-(M 2 K + k-p x + k-p 2 -2k-p 3 

+2p 1 -P2-P1-P3-P2- Pz)B Q (p 3 - k, M w , M n ) 

-{M 2 K -Ml-4k-p 3 )(M 2 K + k- Pl + k-p 2 -2k- p 3 

+2 Pl -P2-P1-P3-P2- P3)C (-p 3 , -k, 0, M w , M w )} . (24) 



Figure 5 
We have 



3 • 



^(8) _ <3 £(27) 

2e 2 



^ = -Z{^ (25) 



^- (Ml + fc • pi + k ■ p 2 - 2k ■ p 3 + 2pi • p 2 - pi ■ p 3 - p 2 ■ p-i) x 
{-2M 2 K (4n)- 2 [l + \n(m 2 ) - \n(M n M K )] 
+A (M 7T )-M 2 w B (0,M n ,M K ) 
+M 2 B (-k, 0, M K ) - 2M 2 K B 1 (-k, 0, M K ) 
-2M 2 K B x (-p 3 , 0, M w ) - (2M 2 K + M 2 )B (-k, 0, M w ) 

+A^(3A4 + A^)C (-A;,-fe,0,M ff ,M^)} , (26) 

2e 2 
zf' ] = —M 2 K (Ml + k-p l + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 -p 3 -p 2 - p 3 ) x 

{-2(4tt)- 2 [1 + ln(m 2 ) - ln(M ff M*)] 

-2S 1 (-p 3 , 0, M ff ) - A,(0, M ff , M*) + B (-k, 0, M ff ) 

+5 (-A;, 0, M K ) - 2B,(-k, 0, M x ) 

+ (3M 2 K + Ml)C (-k,-k,0,M w ,M K )} , (27) 

Z{ 8) = -ZP = | Z{ 27) (28) 

2e 2 
= — — M 2 K (Ml + k-p 1 + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 -p 3 -p 2 -p 3 ) x 

[-MM*) + 2(M 2 K + M 2 n )B (-k, 0, il^)] . (29) 



Figure 6 
We have 



r/ 8) = -t/ 27) (so) 

= -^ {-AoiM^) + A (M K ) - M*B (-k,0,M K ) 
+ (M| + 6M 2 + k- Pl + k-p 2 -2k-p 3 
+2p x -p 2 -Pi-P3-P2- P3)B (-p 3 , 0, M w ) 
-{M 2 K + 5M^ + k- Pl + k-p 2 -2k-p 3 

+2pi -P2-P1-P3-P2- P3)B (-P3, 0, M X ) 

+M 2 K B (p 3 - k, M K , M K ) - M 2 K B (0, M n , M K ) 

-{2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2p 1 -p 2 -pi-p3-P2- P3)B {p3 - k, M v , M K ) 

+{M 2 K + 3M 2 + k- Pl + k-p 2 -2k-p 3 

+2p 1 -p 2 -pi-p3-P2- Ps)B (k - p 3 , M K , M K ) 

+ {M 2 K + Ak ■ p 3 )(2M 2 K + 2M 2 + k- Pl + k-p 2 -2k-p 3 

+2p 1 -p 2 -pi-p 3 -p 2 - p 3 )Co(-P3, ~k, m 7 , M w , M K ) 

+M 2 K (M 2 K + 7 Ml + k- Pl + k-p 2 -2k-p 3 

+2p x -p 2 -pi-p 3 -p 2 - P3)C (-P3, -P3, 0, M v , M K ) 

-(M 2 K -Ml + Ak- p 3 )(2M| + 3M^ + k ■ Pl + k ■ p 2 - 2k ■ p 3 

+2pi -p 2 -pi-p 3 -p 2 - P3)C (-P3, -k, 0, M K , M K ) 

-M 2 K (2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2 Pl -p 2 -pi-p 3 -P2- P3)C {0,p 3 - k, M n , M K , M K ) 

+M 2 K (M 2 K -Ml + Ak- p 3 ){2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p 3 -p 2 - p 3 )D (-p 3 , -k, -p 3 , m 7 , M v , M K , M K )} , (31) 

Tl 8,) = ^-M 2 K {-B (-k,0,M K )-B (-p 3 ,0,M n ) 
o 

+B (p 3 - k, M K , M K ) + 2B {-p 3 , 0, M K ) - B (0, M w , M K ) 

+ (2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p 3 -p2- P3)C (-P3, -k, m 7 , M n , M K ) 
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+{M 2 K + 7M 2 + k-p x + k-p 2 -2k-p 3 

+2p x -P2-P1-P3-P2- P3)C (-P3, -p 3 , 0, M v , M K ) 

-(2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2 Pl -P2-Pi-P3-P2- P3)C {0,p 3 - k, M n , M K , M K ) 

-(M 2 K -Ml + Ak- p 3 )C (-p 3 , ~k, 0, M K , M K ) 

+ (M 2 K -Ml+Ak- p 3 )(2M 2 K + 2Ml + k- Pl + k-p 2 -2k-p 3 

+2p 1 -P2-Pi-P3-P2- P3)D (-P3, -k, -p 3 , rrij, M w , M K , M K )} . (32) 



Figure 7 
We have 



M 2 ms>) _ 3 
Ml A ~ 2 



*(8) _ M n MS') _ * ~(27) 



zr = -T^zr = -:zr m 



e 



Ml {-2A (M K ) 
o 

+ (5M| + 3Ml + k- Pl + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p3-P2- p 3 )B (-k, 0, M K ) 

+ {3M 2 K + 5M 2 + k- Pl + k-p 2 -2k-p 3 

+2pi -p 2 -p\-p 3 -p 2 - P3)B (-P3, 0, M K ) 

-(Ml + 3Ml + k- Pl + k-p 2 -2k-p 3 

+2pi -p 2 -pi-p 3 -p 2 - p 3 )B (p 3 - k, M K , M K ) 

+ {M 2 K -Ml + Ak- p 3 ){M 2 K + 3Ml + k- Pl + k-p 2 -2k-p 3 

+2 Pl -P2-P1-P3-P2- P3)C (-p 3 , -k, 0, M K , M K )} . (34) 



Figure 8 
We have 

Zf ] = \zf 7) (35) 

e 2 
= — (2M 2 K + 2Ml + k-p 1 + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 -p 3 -p 2 -p 3 )x 



6 

{-2Ml(Any 2 [I + ln(m 2 ) - \n(M n M K )] 
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+A (M K ) - 2M*B (-p 3 , 0, M K ) - M 2 K B (0, M n , M K ) 
-2M 2 B 1 (-p 3: 0,M 7r )-2M 2 B 1 (-k : 0,M K ) 
+4M 2 M 2 K C (-p 3 ,-p 3 ,0,M n ,M K )} , (36) 

2 

zf' ] = --M 2 K (2M 2 K + 2M 2 + k-p 1 + k-p 2 -2k-p 3 + 2p 1 -p 2 -p l -p 3 -p 2 -p 3 )x 
o 

{-2(4tt)- 2 [1 + ln(m 2 ) - HM W M K )] 

+2B (-p 3 , 0, M K ) - B (0, M w , M K ) - 2B 1 (-p 3 , 0, M w ) 

-2B x {-k^,M K ) + AMlC Q (-p 3 ,-p 3 ^,M^M K )} , (37) 

Zr = -^ZP = 3 -Zr (38) 

e 2 
= --Ml (2M 2 K + 2Ml + k-p 1 + k-p 2 -2k-p 3 + 2p l -p 2 -p l -p 3 -p 2 - p 3 ) x 
o 

[-A (M K ) + 2(M| + Ml)B (-p 3 , 0, M*)] . (39) 



Figure 9 



We have 



M, 



2 



Zf> = -^Z?> = -Zr> (40) 



C-: 



L K 

J M 2 {-2A (M X ) 
-2(ft ■Pi + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 -p 3 -p 2 - p 3 ) x 
(4tt)- 2 [1 + ln(m 2 ) - ln(M ff Mtf)] 
+ (8M 2 + k ■ pi + fc • p 2 - 10fc • p 3 
+2pi • p 2 + 3pi • p 3 + 3p 2 ■ p 3 )B (-p 3 , 0, Af ff ) 

-2(k ■pi + k-p 2 -2k-p 3 + 2p 1 -p 2 -p 1 - p 3 - p 2 ■ p 3 )Bi(-p 3 , 0, M n ) 
+(8M| -3k-p 1 -3k-p 2 -2k-p 3 
+2pi -p 2 -pi-p 3 -p 2 - p 3 )B (-k, 0, M x ) 

-2(k-pi + k ■ p 2 - 2k ■ p 3 + 2p ± ■ p 2 - p ± ■ p 3 - p 2 -p 3 )Bi(-k,0,M K ) 
-(2M 2 + k- Pl + k-p 2 -Ak-p 3 
+2pi -p 2 -pi-p 3 -p 2 - p 3 )B (p 3 - k, M w , M K ) 
+2(-k ■ pi - k ■ p 2 + pi • p 3 + p 2 ■ p 3 )B 1 (p 3 - k, M n , M K ) 
+4(fc • p 3 ) (k ■ p ± + k ■ p 2 - 2k ■ p 3 
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+2pi -p 2 -pi-p3-P2- P3,)Co(-p 3 , -k, m 7 , M n , M K ) 

+8(k ■ p 3 ){ Pl -P3+P2- Pz)Ci{-ps, -k, 0, M w , M K ) 

+8(k .p 3 )(k- Pl + k- p 2 )C 2 (-p 3 , -k, 0, M n , M K )} . (41) 

Figure 10 

We have 

r/ 8) = ~ T/ 27) = -(M| - M 2 )- 1 ^ 8 > = -| {Ml - M 2 )" 1 ^ = . (42) 

This result is somewhat expected since the neutral pion is electrically neutral. 

Our amplitude is ultraviolet divergent. So is the amplitude calculated in Ref. [2|. Only 
the sum of these is ultraviolet finite. Our amplitude is also infrared divergent. We shall 
show in the following section that the infrared divergence cancels the one coming from real 
soft photon emission. The cancelation occurs at the level of differential decay rates. Since 
we are interested in evaluating the amplitude, we have to subtract infrared divergence. The 
subtraction scheme is not unique. We will adopt a minimal subtraction scheme consisting 
on dropping out only ln(m 2 ) terms from the expression of the amplitude. 

IV. INFRARED DIVERGENCE 

The infrared divergence of the amplitude reads: 

2 



e 



Mm -- -i-i Mo [l~(k- pM-Vz, ~k, M ff , M K )} ln(m 2 ) . (43) 



The divergence is canceled by the one generated from the emission of a real soft photon. 
The energy q of the latter being raised by the detector resolution u>. The cancelation takes 
place as usually at the differential decay rate level. Let .M 7 denotes the decay amplitude 
corresponding to the emission of one real photon. We shall focus only on the infrared 
divergent part of M. 1 . The Feynman diagrams one needs are drawn in Fig. 11. The result 
is standard: 

Mi = ^ eMo (t^ - a^) . (44) 

\k-q p-i-q J 
One obtains after squaring the amplitude and summing over photon polarizations 



M n | 2 = -e 2 |M) r 



M\ 2(k-p 3 ) Ml 



+ 



(k-q) 2 (k-q)(p 3 -q) (p 3 ■ q) 2 _ 
12 



(45) 



We shall now explicitly show the infrared divergence cancelation. The infrared divergence 
coming from virtual photons is deduced from Eqs. (J3J) and (J4*3*|) 



e 



2 



dT IR = -— dT [l-(k- Pa)r(-P3, -k, M w , M K )\ ln(m 2 ) . (46) 

From Eq. (|45|). the differential decay rate corresponding to the emission of one real photon 
reads in the soft photon approximation 

(47) 



(48) 



d^ = e 2 dr r , d \i 


\ Ml 


2(k-p 3 ) Ml 
(k-q)(p a -q) (p 3 ■ q) 2 . 


_(k-q) 2 


One can then easily check that 




dT IR 


+ dT] R 


= 0. 



V. RESULTS 

In this section we shall evaluate isospin breaking correction due to explicit virtual pho- 
tons. This correction depends only on S3 thanks to the on-shell relation (J3J). It is natural 
then to study the variation of the correction with respect to S3 in the physical region defined 
by Eq. (JBJ). In order to evaluate the size of the correction it is convenient to compare it to the 
tree-level value or to the one-loop level correction in the absence of isospin breaking. Since 
the latter depends on two kinematical variables it is easier (but not less instructive) to com- 
pare it to the former. This is done in Fig. (12) where we plotted the Born amplitude squared 
I M. | 2 (dashed curve) and the one-loop level amplitude squared | M.q | 2 +2 3? (M.qM.\) 
(plain curve). We recall that M.\ stands for the explicit virtual photon correction. In order 
to draw Fig. (12) we used the following numerical input: 

e 2 = (4vr)/(137.036) , M w = 0.134977 GeV, M K = 0.495042GeV, 

C = -1.07 x 10- 6 GeV~ 2 , G 8 = 5.45, G 27 = 0.392. (49) 

Note also that the ultraviolet divergent terms (A terms) have been dropped out and the 
scale /1 was taken to be 0.770 GeV. Let us comment the content of Fig. (12). To this end, 
we shall write the amplitude squared as 

I .M | 2 = I .Mo I 2 [l + ^s)] • (50) 

The ratio 5 is given in Tab. |IJ for different values of S3. As can be seen from Tab. (jlj), the 
explicit virtual photon correction is about 5 to 6% the value of the Born amplitude. 
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VI. CONCLUSION 

In this work we calculated isospin breaking corrections to the process K + — > 7r + 7r°7r° 
due to virtual soft photons at one-loop level and in the framework of chiral perturbation 
theory. The corrections generated by m u — m^ and virtual hard photons have been calculated 
in Ref. 2J. The corrections are individually ultraviolet divergent but jointly ultraviolet 
finite. They are also infrared divergent. We showed that this divergence is canceled at the 
differential decay level by the one coming from real soft photons. We follow Ref. |2j and 
denote by A 00+ (0, 0) the decay amplitude for the process in question evaluated at the center 
of Dalitz plot. Then the amplitude squared reads: 

I A 00+ (0,0) | 2 = 2.49 x i(T 13 + 6.84 x KT 13 



o( p 2) e>(p 4 ) 

+ 7.7 x l(T 14 + -1.3x KT 14 , . (51) 

is. br. I is. br. II 

In the first line of the equation we reported the value of the amplitude squared at leading 
and next-to- leading chiral orders and in the absence of isospin breaking 5]. The first term 
in the second line represents isospin breaking corrections due to m u — rrid and to hard virtual 
photons J2j. The second term has been calculated in the present work and corresponds to 
the isospin breaking correction due to soft virtual photons. Adding all this together, isospin 
breaking correction to the amplitude squared represents 7% its one-loop value in the absence 
of isospin breaking and at the center of Dalitz plot. 



APPENDIX A: LOOP INTEGRALS 



We use dimensional regularization and adopt the MS subtraction scheme 



A 



32vr 2 



1 - 7 + ln(47r) 



(Al) 



A-n 

where n is space-time dimension and 7 the Euler constant. All the technical material 
necessary for the calculation of one- loop integrals is given in the appendix of Ref. jq|. 
It is convenient to take the following notations: 

/ 4M 2 /k s 

ar = Jl--f., (A2) 
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A integrals 

The one-point function reads: 

A (m) = m 2 
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1 . i m* 
In 
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16tt 2 



where \x an arbitrary scale with mass dimension. 



B integrals 
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We need the following two-point functions: 
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r integrals 

The definition and expression of r integrals can be found in Ref. 
expression 
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C integrals 



The three-point functions needed for our purposes are: 
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C 2 {-p 3 ,-k,0,M n ,M K ) = 

\-\ Sz ,MlM 2 K ) [-2M 2 w B {-p 3 ,0,M n ) 
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D integrals 

We have one complex and one real four-point function in the physical region. The complex 
one reads: 
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The real one is given by: 
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r 


1 


1.12 


1.24 


V(4M2) 


1.48 


1.60 


(M K - Mrf/{AMl) 


10 2 <5 


-5.78 


-5.59 


-5.43 


-5.26 


-5.13 


-4.99 


-5.44 



TABLE I: The size of explicit virtual photon correction compared to the Born amplitude. We used 
the notations r and sq for the ratio s/(4M^) and the Dalitz plot center (S1 + S2+S3 = Mj^+3M^)/3, 
respectively. 
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Figure 1: Born diagrams. The strangeness-changing vertex is represented by 
a full squarre. 
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Figure 2: The various topologies with a virtual photon (wavy line) attached 
to the K'K'K'K strangeness-changing vertex. 
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Figure 3: Box diagrams with a 7T7t scattering vertex. 
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Figure 4: Diagrams with a pion pole and a photon attached to a 7T7t scattering 
vertex. 



x^h ^_* n+ k+ » , r^ .. » 



Jf + > J ■ ' A > 7T 





7T° ^° 7T° 7T° 

(a) (b) 

^ + > r^^ „ 7T+ X+ > J— 2 > 





7T° ^° 7T° ^° 

(c) (d) 

#+ > r^ ■ ., > 7T+ X+ > ■ _ P^ > 7T+ 





(e) (f) 



Figure 5: Diagrams with a pion pole, a nn scattering vertex, and a photon 
attached to the Kit strangeness-changing vertex. 
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Figure 6: Box diagrams with a .fOr scattering vertex. 
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Figure 7: Diagrams with a kaon pole and a photon attached to a i^7r scat- 
tering vertex. 
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Figure 8: Diagrams with a kaon pole, a i^7r scattering vertex and a photon 
attached to a Kir strangeness-changing vertex. 
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Figure 9: Diagrams with a neutral kaon pole. 
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Figure 10: Diagrams with electromagnetic pion form factor. 
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Figure 11: Emission of one real soft photon diagrams. Are shown only dia- 
grams generating infrared divergence. The filled circle represents tree vertices 
from Fig. 1. 
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Figure 12: The amplitude squared as function of S3. The dashed curve 
represents Born amplitude squared | M.q | 2 , the plain one includes explicit 
virtual photon correction | Mq | 2 +2$t (M.qMi) . 



